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Abstract. Entanglement dynamics of two nonintcracting qubits, locally affected by 
random telegraph noise at pure dephasing, exhibits revivals. These revivals are not due 
to the action of any nonlocal operation, thus their occurrence may appear paradoxical 
since entanglement is by definition a nonlocal resource. We show that a simple 
explanation of this phenomenon may be provided by using the (recently introduced) 
concept of hidden entanglement, which signals the presence of entanglement that may 
be recovered with the only help of local operations. 
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1. Introduction 

Entanglement is one of the most peculiar features of quantum mechanics and it 
also plays the role of a fundamental resource in many applications of quantum 
information [TJ [2j [3j H] . On the other hand, entangled systems unavoidably interact with 
their environments causing decoherence and a loss of entanglement. Since entanglement 
is by definition a nonlocal resource, one expects that any attempt to restore it must 
involve the use of nonlocal operations. 

We consider physical situations where two subsystems, for example two qubits, are 
prepared in an entangled state and subsequently decoupled [SI IB]- Due to the interaction 
with their local environment, entanglement dynamics may exhibit a non- monotonous 
behaviour, with the occurrence of revivals alternating to dark periods [7| |8| l9l IT0 | ITT | [T2] . 
In some cases, this phenomenon is due to the fact that entanglement is transferred to 
quantum environments, and then back-transferred to the system [71 [8]. In other cases, 
the environment can be modelled as a classical system [9j [HJ [13] and no entanglement 
between the system and the environment is established at any time. In the latter cases, 
the occurrence of entanglement revivals may appear paradoxical, since the effect of 
the noise is analogous to a local operation on a subsystem. A first interpretation of 
this phenomenon has been given in terms of correlations present in a classical-quantum 
state of environments and qubits [TT]. In a recent work [13] . we have proposed to solve 
the apparent paradox by introducing the concept of hidden entanglement (HE), which 
measures the amount of entanglement that may be recovered without the help of any 
nonlocal operation. The definition of HE is based on the quantum trajectory description 
of the system dynamics, that allows to point out the presence of entanglement in the 
system even if the density operator formalism does not reveal it: this entanglement is 
thus not accessible (hidden) due to the lack of classical information [14] . 

Relevant examples of situations where the environment can be modelled as a 
classical system may be found in solid-state implementations of qubits. For example, in 
superconducting nanocircuits one of the most relevant sources of decoherence [T5l[T6l [T7] 
are fluctuating background charges localized in the insulating materials surrounding 
superconducting islands [H]. Each impurity produces a bistable fluctuation of the 
island polarization. The collective effect of an ensemble of these random telegraph 
processes, with a proper distribution of switching rates, gives rise to l//-noise [20] 
routinely observed nanodevices [HI [18], fT9] . 

In this paper we exploit the concept of hidden entanglement to explain the 
occurrence of entanglement revivals in a simple system. In particular, we consider 
two noninteracting qubits, one of them affected by a random telegraph noise at pure 
dephasing [211 1221 [23] . The paper is organized as follows. In section [21 we introduce the 
Hamiltonian model. In section [3] we discuss the entanglement dynamics, showing that 
the revivals of entanglement are due to the presence of hidden entanglement. In section 
HJwe summarize obtained results and present some final comments. 
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2. Model 



We consider two noninteracting qubits A and B, initially prepared in a pure maximally 
entangled state p(0) = \<p)(<p\, evolving according to the Hamiltonian {h = 1) 
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where a ZA = a z ® 1, o~ ZB = 1 ® °~z an d STi represents a random telegraph (RT) process 
G {0,w} [23] acting on qubit A. The RT process induces a random switching of 
qubit A frequency between Qa/(2tt) and (Qa + u )/(2tt), with an overall switching rate 
7 (without loss of generality, we assume v > 0). We consider a symmetric RT process 
where the transition rates between the two states are equal, that is 7o->„ = Jv^-o = 7/2. 

Our first aim is to find the system density matrix at any time p(t). The two qubits 
independently evolve under the Hamiltonian ([T]): qubit B freely evolves whereas qubit 
A displays a pure dephasing dynamics due to the effect of the stochastic process 
The dynamics of single qubit subject to RT noise at pure dephasing has been solved 
. A possible way to obtain p(t) is to solve a stochastic Schrodinger equation 



in 



which gives the following formal expression for p(i) 



p(t)= V[t(t)]P[ffl\ P&) 



(2) 



where p ? (t) = \if^(t))(ip^(t)\ with \(pt(t)) = ea Jo inot \(p), and the probability of 

the realization can be written as 

P[£(t)] = lim 

m— >oo 

where r] m+ i is a (m + 1) joint probability for the sampled at regular intervals 
At = (t — to)/m, tk = to + kAt, (,k = £{tk) (k — 0, . . . , m) [25]. Since the qubits evolve 
independently, the above procedure leads to a simple form depending on the single qubit 
coherences. In the computational basis {|00), |01), |10), |11)}, where \ij) = \i)<8>\j), with 
&z\i) = (— an d i £ {0, 1}, and assuming an initial Bell state \ip) = ^(|00) + |11)), 
we obtain 

/ 1 q (t)e l ^ A+nB ^ \ 
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where the coherence decay factor q{t) reads [2D 



P(t) 
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with a = \J\ — g 2 , A — |(1 + -) and g — v/j. In the following we shall exploit pit) 
given by Eq.© to analyze the two-qubit entanglement dynamics. 
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3. Entanglement dynamics 

To quantify the degree of entanglement of the system state p(t) we use the entanglement 
of formation Ef [26\ that can be readily calculated by the formula [27] 

where C(p(t)) is the concurrence and h(x) = —x log 2 x — (1 — x) log 2 x(l — x). For the 
state p{t) of Eq. (j4j) we obtain C(p(t)) = \q(t)\, where q(t) is given in Eq. ([5]). It is 
worth to notice that the evolved state p(t) belongs to the Hilbert space spanned by the 
Bell states \4>±) = -^(JOO) ± |11)), therefore the entanglement of formation equals the 
entanglement cost [2B]. In the strong coupling regime, g — v/j > 1 [211 123] entanglement 
revivals occur during the system dynamics [9j [T3] . 

3.1. Dephasing under a "static" RT process 

To understand the nature of this phenomenon we initally consider the limiting case of 
an extraordinarily slow RT process, 7 — > (g — > 00). This regime physically describes 
situations where the stochastic process is slow enough to be considered static during the 
system time evolution lasting t, i. e. we assume I/7 3> t [To]. The evolution expressed 
by Eq.(j4]) describes an average resulting from the collection of several time evolutions 
each lasting t. The average includes the possibility that the RT process takes any of the 
two values £ = or £ = v at time t = with equal probability. By a straightforward 
calculation we find that the concurrence in this case is given by C(p(t)) = | cos(vt/2)\. 
Under these conditions we do not find any entanglement decay, rather the concurrence 
is equal to one at times t n = 2mr/v and it vanishes at times t n = (2n + l)ir/v , where n is 
a non-negative integer. The entanglement revivals (see the top solid curve in Fig. [TJ) are 
not due to periodic entanglement death and rebirth by nonlocal operations. Indeed, the 
Hamiltonian evolution described by Eq. ([T]) only includes local operations. Since local 
operations cannot increase entanglement [3111], its increase during the intervals ]£„,£„] 
must be attributed to the manifestation of quantum correlations that were already 
present at times t n , but were hidden, in the sense that the density operator formalism 
does not capture them. 

These correlations are evident in the quantum trajectory description of the system 
dynamics [29]. The system evolution in fact results from averaging on only two possible 
quantum trajectories. The first trajectory corresponds to = and the Bloch vector 
of qubit A rotates around its z-axis with frequency Qa/{^)- The second trajectory 
orresponds to — v and the Bloch vector of qubit A rotates around the z-axis with a 
different frequency (Qa + ' i; )/(2 7r )- The Bloch vector of qubit B instead rotates in both 
cases around its z-axis with frequency Q b /(2ti). Thus, during the first trajectory the 
system state evolves, up to an irrelevant global phase factor, as 

|^ (t)) = i=(|00)+e^ + ^ t |ll)), (7) 
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Figure 1. (color online) Entanglement of formation Ef(p(t)) as a function of the 
dimensionless time vt for the \tp) = -i=(|00) + |11)). From top to bottom, the solid 
curves correspond to g = oo, 200, 50, 10 and 5, the dashed curves represent /(e~ 7 */ 2 ) 
for the same values of 7 = v/g. 



while during the second trajectory the system evolves, apart from an irrelevant global 
phase factor, as 



|V«(*)> 



W|00) + e- ivt e- i ^ A+slB)t \ll)y 



s/2 



(8) 



The two quantum trajectories only differ by the fact that the basis states 1 00) and |11) 
they acquire the additional relative phase ivt in the quantum superpositions of Eqs. flTj) 
and flS]). Since the two quantum trajectory occur with equal probability, the system's 
state is described by the quantum ensemble 

A = {(p ,\Mt))),(pv,\Mt)))}, (9) 

where Po = p v = \- The entanglement associated to the quantum ensemble A can be 
suitably identified by its average entanglement given by [261 USB EH E21 [33] 

E av (A,t) = E PiE(\ipi(t))) = l, (10) 

i£{0,v} 

since both states \(fo(t)) and \ip v (t)) are maximally entangled at any time (E is the 
entropy of entanglement [31 Hj). 

The hidden entanglement [H] of the ensemble A is defined as the difference between 
the average entanglement of Eq. ([TO]) and the entanglement of the corresponding density 
operator p(t) = Y*Pi\ < Pi( t )){ < Pi( t )\ : 

E h (A,t) = E av (A,t)-E f (p(t)). (11) 

The meaning of hidden entanglement is simple: it measures the entanglement that 
cannot be exploited as a resource due to the lack of classical knowledge about which 
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state in the ensemble A we are dealing with. Once this classical information is provided, 
the entanglement can be recovered. 

Coming back to the interpretation of entanglement revivals, the ensemble 
description (average entanglement E av (A,t) = 1) tells us that at times t n the system 
is always in a maximally entangled state {\<fo(t n )) or \<p v {t n ))) but the lack of classical 
knowledge about which of the two states in the ensemble A we are dealing with prevents 
us from distilling any entanglement: in fact, entanglement is hidden being Ef(p(t n )) = 
and Eh(A(t n )) = 1. At times t n this lack of knowledge is irrelevant since the random 
relative phase becomes meaningless at vt n = 2nn: all the initial entanglement is 
recovered, Ef(p(t n )) = 1 and Eh(A(t n )) = 0. 

3.2. Dephasing due to a RT process: dynamic case 

We now investigate the case when the RT process evolves during the systen evolution 
time, i. e. we consider the regime I/7 > t. This situation is illustrated in Fig. [T] where 
we observe that the amplitude of revivals decreases as 7 increases (g decreases). Aslo in 
this case there is hidden entanglement. The possible quantum trajectories the system 
undergoes are now infinite. The system state is described by the quantum ensemble 
A(t) = {P[C(t}], \(P((t))} and the average entanglement of A is calculated by solving the 
path-integral 

E av (A,t) = J V[ffl]P[ffl]E(\<P&)))- (12) 

Once again, we obtain E av (A,t) = 1 since during each trajectory the state remains 
in a maximally entangled state at any time. On the other hand, the entanglement 
of formation assumes lower values with respect to the static noise case, 7—7-0. In 
particular, the amplitude of revivals does not reach anymore the initial maximum value. 
This is due to the fact that, in general, the action of the RT process during the time 
evolution makes the coherences of |^(t))(v 5 5(^)| (i n the basis {|00), |11)}) no longer in 
phase at times t n . In the time interval ]0, t n ] one or more transitions can occur between 
the two RT states, such that we can have a random extra phase at the times t n given 
by 

■&{t n ) = / dt'iit') - 2vm, (13) 
Jo 

where n is a non-negative integer. This unknown phase difference is responsible for 
the decay of the absolute values of coherences |g(t n )| in the evolved two-qubit state 
p(t n ) of Eq. (HI): if we knew the phase difference "d(t n ) for each state \<p^(t)), we would 
be able to restore the coherence absolute value to 1, and therefore recover all the 
initial entanglement, by simply applying the unitary local operation e 1 2 °*a. For 
completeness we point out, that the relative maxima of the entanglement of formation 
occur at t* = t n j — jp, as one can derive from equation ([5]). 

Notice that the amount of the decay of the amplitude of entanglement revivals is 
monotonously related to the normalized autocorrelation function of the symmetric RT 
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process R(r) = (£(t)£(0))/(£ 2 (0)) = e~ 7 * [21]. Indeed, as we have already mentioned, 
the reduction of the amplitude revivals is related to the transitions of the RT occurred 
in [0, t n ], whose mean number is yt n /2. From a quantitative point of view, for g > 1 
the coherences decay factor Eq. (jSJ) can be approximated as \q(t)\ ~ e~ 7< / 2 [cos(-ut/2) + 
l/gsm(vt/2)], so that C{p{t n )) ~ e" 7 *™/ 2 and E f (p(t n /2)) = f(-~/t n /2), with / defined 
in Eq. ([6]). This clearly shows that the decay of the entanglement amplitude revivals is 
due to the decrease of the RT correlations, or in other words, to the memory loss of the 
stochastic process £(t) itself. 

4. Conclusions 

In this paper we have exploited the concept of hidden entanglement to interpret the 
occurrence of entanglement revivals in a particular system where back-action from 
the environment is absent. Namely, we have considered a system composed of two 
noninteracting qubits where one qubit is subject to random telegraph noise at pure 
dephasing. During the dynamics, entanglement vanishes and revives always "remaining" 
inside the system, as it is signalled by the average entanglement, E av (t) = 1 at any 
time. At certain times t n this entanglement is completely hidden, in the sense that 
the entanglement of formation Ef(t n ) = while the hidden entanglement Eh(t n ) = 
E av {tn) — Ef(t n ) = 1. For this reason, the two-qubit entanglement can be simply 
recovered at subsequent times without the help of any nonlocal operation: in the 
considered case in fact the Hamiltonian only involves local operations. 

Finally we remark that the concept of hidden entanglement can be of practical 
relevance in solid-state devices, where dominant noise sources tipically have large 
amplitude components at low frequencies. In these systems entanglement revivals may 
also be induced by applying local pulses to the qubits [331 El El] • 
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